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(RADIUS ARTERY 
PULSEWAVE) 



(54) LIVING BODY CONDITION MEASURING APPARATUS 

(57) An apparatus for uninvasively evaluating hemo- 
dynamic parameters, particularly a sphygmographic 
analyzer for evaluating vessel compliance and flow 
resistance in the central and peripheral portions of the 
artery system separately. A microcomputer (4) detects 
the radial artery wavef rom of an object through a sphyg- 
mograph (1), and fetches a stroke output of the object 
from an output meter. Next, the value of each of the five 
elements that constitute an electric circuit simulating the 
artery system ranging from the central portion to the 
peripheral portion of the living body is adjusted on the 
basis of the stroke output so that when an electric signal 
corresponding to a pressure wave of the aortic origin of 
the object is applied to the electric circuit, its response 
waveform corresponds to the radial artery waveform, 
and the resulting value of each element is outputted as 
the hemodynamic parameter. Further, the pulse wave- 
form, the maximum blood pressure value and the mini- 
mum blood pressure value at the aortic origin are 
calculated from the value of each element, and the cal- 
culation result is outputted to an output device (6). 
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Description 

1. Field of the Invention 



The present invention relates to an optimal device for measuring physiological state which is used to measure con 
£££ J« ^ ^ SPeC * iCa " y - 816 Present invention relates to « <**• for a^zin^X «ves T^e 

H C ^ t0ry T 6 " 1 in 1,16 hUman b0dy by evaluatin 9 the central blood Pressure and tneSn^Sce and 
resistance in the blood vessels at the central and peripheral portions of the circulatory system. comp,,ance and 

2. Background Art 

HoJ^£ r nT!t ana ^ rate m0St common, y used for th « diagnosis of the condition of the circulatory system 
However. ,n order to carry out a more detailed diagnosis, it is necessary to measure the so-called circulatory ^me- 
ters of compliance and viscous resistance in the blood vessels circulatory parame 
The pressure waveform and blood flow volume of the proximal portion of the aorta and at the site of insertion of a 
catheter into a peripheral artery need to be measured in order to measure these circulatory parameters or an indir^ 

^STS ^ be 3PP,iCable HOWever ' the former method is an'inva^e one. and emp^ 
Jf « ' "5 6 thC ' atter meth0d ' althOL,gh permittin 9 ^n-^vasive measurement of blood flow wiSS 

A^T^T U ' reS fra ' n6d technicians trai ™9 and, moreover, necessitates a large and expens-ve^cT 
fournt?' ^ ^Si '""^ d,scovered a "ethod '°r approximating the parameters accor*ng to a lumpS 
ZlZ S by measuring only the pulse waveform at the radial artery and the stroke volume T^SeaSne 

present inventors proposed a pulsewave analysis device capable of carrying out an evaluation of the crrculatorvolrJm 

Sy^ 

t^n^^V* 6 afor « mentioned method does not take into account the compliance at the peripheral and central por- 

JSSSTE!-?!^ Hei 4 ' 276234 ' diSCl0ses a sphygmomanometer at the periphery of^paT- 

tients body. Namely, as shown in Fig. 29, cuff 1 10 is wrapped around the upper arm of a patient and a band 138 is 

oSZZ J 6 Pat t!, S W T1 14 °- PulS6WaVe S6nS0r 134 is putted °" * e radiaT arte? ^ oHhe %LTand the 

, ,r J^T'' * bl0Od Pr6SSUre 3t the P eri P her V and ^nter of the arterial system in the human body are actually meas- 

mofe the dZSof SS^^^S^ " in the ^o^ P rl^lZt 

more, the degree of this difference vanes depending on the shape of the pulsewave of the perioherv of the arterial 
system. Figs. 22 through 24 are provided to explain this variation in blood pressure depend^ on £ ^^1^ sSe 
The pressure waveform and systotic/diastolic pressure of the aorta, which is central portion of ^^2^^ Tnd 
the pressure waveform and systolic/diastolic pressure at the radial artery are shown in these figures ^ 

the dlshS nn^nl uhf ?r ° f PU ' Se WaVef ° rm ' WhSrein the SyStolic pressure obtained at *e aorta is indicated by 
the dashed line and the systolic pressure obtained from the radial artery is indicated by the solid line Althouah the 

^ ^ T^l ^ " S,i9htiy hi9her ' *- b ' 00d preSSure ^ ^ -^most XiSa 
£L th6 h ? Se ° f tne u second ^e of P"'se waveform shown in Fig. 23, however, the difference between the systotic 

ZTZl t , at , the a0rta ^ 3t radia ' art6ry iS 14 9 mmH 9' a considerable difference thaHLrv^ in me 
Ra 24 H«! P k ?T Sh ° Wn " R9 ' 22 " FurthSr ' in 1,16 C3Se 0f the third * pe of P ulse reform shown in 
W 4 'n . T 6lWe ! n 1,16 SySt0 ' iC preSSUre iS 9reater ' 26 - 1 mmH 9- Mo ^over, in contrast to the TyLTand 

ISrlit TthT • ' n *! 0356 °' TyPS '" PU,S6 Waveform ' 106 pressure wavefo ™ of *e aorta is higher in its 

not S^h . ? h Pr6 f ?u ef ° rm ° btained 3t th6 radial arter * However ' the diast ° lj c Pressure at the radiafartej do 
not depend on the shape of the pulsewave, but are approximately same 

t^UTh ' P r S lT.!! 0r T iS ° bSerVed in 3 n ° rmal healthy person Tne waveform is relaxed and loose, and is charac- 
terized by a fixed rhythm by regularity, with little disruption. On the other hand. Type II pulse waveform shows a shTm 
nse. The dicrotic notch is deep, and the subsequent peak of the ejection wave is higher the usual TypeC^e wav£ 

TrhSthS flT P T Ure r T inin9 e,eVat6d f ° r ,0n96r P6ri0d ,han in Types 1 and " wSefo 
elev^ti ^Jth k, J, 9 • 11 ' S POSS ' ble <0r the P6ripheral bl00d pressure of the radial artery or upper arm to be 
elevated, while the blood pressure of the proximal portion of the artery, i.e., of aorta is low. Further the opposite srtua 
tion is also possibte, namely, the blood pressure of the periphery is lower than that of the central po L^e artena. 
system ,s high. These differences are depending on the shape of the pulse waveform 
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For example, when an anti hypertensive agent is administered to a patient, it is possible that the blood pressure at 
the central portion of the arterial system is not actually reduced, even if there is a drop in the blood pressure of the 
periphery. It indicates that it is not reliable to ascertain drug effect based only on the peripheral blood pressure. Con- 
versely, even if no change is observed in the blood pressure of the peripheral portion of the arterial system, the actual 
5 load on the heart may in fact have been reduced if pulse waveform changes to type II shape. As a result, it may be con- 
cluded that the conventional method to measure peripheral blood pressure is not reliable to evaluate the true cardiac 
load by the central blood pressure. 

DISCLOSURE OF THE INVENTION 

10 

The present invention was conceived in consideration of the above circumstances, and has as its first objective the 
provision of a device capable of evaluating circulatory state parameters in a non-invasive method, and more particularly, 
to the provision of a device for analyzing pulse waves that is capable of evaluating compliance and resistance in blood 
vessels at the center and periphery of the arterial system separately. 
15 Further, the present invention's second objective in the provision of a sphygmomanometer capable of obtaining the 

blood pressure at the center of the arterial system from pulse waveforms measured at the periphery of the arterial sys- 
tem. 

The first standpoint of the present invention is characterized in the provision of a measuring means for measuring 
physiological slate, and an analyzing means for calculating circulatory state parameters, including the viscoelasticity of 
20 the aorta, as circulatory state parameters which show the circulatory state of an organism s arterial system from the 
center to the periphery thereof, based on the physiological state of the organism. 

Thus, since circulatory state parameters, including viscoelasticity at the aorta, in the arterial system of an organism 
are calculated, it is possible to achieve a more precise evaluation since the circulatory state of the arterial system, which 
extends from the center to the periphery of the organism's body, is separated into a center component and a periphery 
25 component. 

In one embodiment of the present invention, circulatory state parameters are determined by approximating the 
aforementioned circulatory state parameters in the organism's arterial system with an electric circuit based on a lumped 
five parameter model. As a result, it is possible to more easily calculate the circulatory state parameters conforming to 
conditions in the human body, as compared to a lumped four parameter model. 

30 In another embodiment of the present invention, the pulse wave at the periphery of the arterial system is employed 
as the physiological state, the pressure waveform at the left cardiac ventricle in the body is provided to a lumped five 
parameter model, and each of the elements making up the lumped five parameter model is determined so as to obtain 
the pulse waveform at that time. As a result, it is possible to determine circulatory state parameters which closely match 
the pulsewave at the periphery which was actually measured from the body. 

35 in another embodiment of the present invention, an electric signal corresponding to the pressure in the left cardiac 
ventricle is approximated by a sinusoidal wave. Thus, it is possible to determine circulatory state parameters which even 
more accurately express current conditions at the center the arterial system in the body. 

The second standpoint of the present invention is characterized in provision of a measuring means for measuring 
physiological state at the periphery of the arterial system, and a blood pressure calculating means for determining cir- 

40 culatory state parameters showing the circulatory state in the arterial system based on physiological state, and calcu- 
lating the pressure waveform at the aorta based on the aforementioned circulatory state parameters. 

Since the present invention calculates the pressure waveform at the aorta from the physiological state at the 
periphery of the arterial system in this way, the blood pressure value at the aorta can be estimated using just conditions 
at the periphery of the arterial system. As a result, there is no concern that a false conclusion will be reached, such as 

45 deciding that an administered pharmacological agent had no effect because no change was observed in blood pressure 
measured at the periphery. Thus, it becomes possible to correctly ascertain the effect of an administered pharmacolog- 
ical agent, making the present invention extremely useful when selecting drugs such as hypertensive agents and the 
like. 

Further, since one embodiment of the present invention provides that parameters including viscoelasticity at the 
so aorta are selected as circulatory state parameters, it is possible to make an evaluation of the circulatory state at the 
center and at the periphery of the arterial system separately. As a result, an even more accurate estimation of blood 
pressure is possible. 

Moreover, another embodiment of the present invention provides that circulatory state parameters are determined 
by approximating the circulatory state of the arterial system using an electric circuit based on a lumped five parameter 
55 model. Therefore, it is possible to obtain the blood pressure value at the aorta which more closely conforms with con- 
ditions in the human body. 

Another embodiment of the present invention employs the pulsewave at the periphery of the arterial system as the 
physiological state, provides the pressure waveform at the left cardiac ventricle to a lumped five parameter model, and 
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determines each element making up the model so as to obtain the pulsewave at this time. As a result it is possible to 

iSZT^-^T? 31 USin9 drcu,ator y state Parameters which closely match the pulsewave at the 

periphery which is actually measured in the body. 

In another embodiment of the present invention, blood pressure at the aorta is estimated by determining circulatory 
state Parameters without detecting stroke volume. Thus, blood pressure measurements can be made without causing 
the test subject any unpleasant discomfort or inconvenience. 

th. errlbod i men ' of the P resent invention, the values of the circulatory state parameters are adjusted so that 

the calculated value of stroke volume which is obtained from the pressure waveform of the aorta is equal to the actual 

me f * S,r ° ke VO ' Ume Which is obtained ,rom * e bod * Thus, it is possible to even more precisely esti- 

mate blood pressure at the aorta. 

In another embodiment of the present invention, the workload on the heart is calculated based on the pressure 
waveform at the aorta. Thus, even when no notable change is observed in the blood pressure value obtained at the 
periphery of the arterial system, it is possible to quantitatively indicate the actual load on the heart. Accordingly, it is pos- 
sible to carry out an even more exact evaluation of a treatment method employing a hypertensive agent, for example. 

BRIEF DESCRIPTION OF THE FIGURES 

Fig. 1 is a block diagram showing the structure of a pulsewave analysis device according to one embodiment of the 
present invention. 

Fig. 2 is a diagram showing the conditions under which measurements employing a pulsewave detector 1 and a 
stroke-volume measurer 2 were conducted in this embodiment. 

huma?body a) * 3 diaQram showin 9 a *>ur parameter model which models the arterial system in the 

humSbod? ^ 3 drCUit dia9ram snowin 9 a lum Ped fiv e parameter model which models the arterial system in the 

th» r«v;l?i I I™??.* 6 pressure waveform at the left cardiac ventricle and the blood pressure waveform at 

Fig. 5 is a diagram showing a waveform which models the blood pressure waveform at the proximal portion of the 

Fig. 6 is a flow chart showing an overview of the operation of the pulsewave analysis device in this embodiment 
Fig. 7 is a ow chart showing the processing operations for parameter calculation in the pulsewave analysis device 
HQ. 8 is a f low chart showing the processing operations for parameter calculation in the pulsewave analysis device 
Hg. 9 is a flow chart showing the processing operations for calculating a and od in the pulsewave analysis device 
Hg. 10 is a flow chart showing the processing operations for calculating © in the pulsewave analysis device 
Hg. 1 1 is a waveform diagram showing an example of the waveform at the radius artery obtained through an aver- 
aging process in the pulsewave analysis device. 

Fig. 12 is a waveform diagram showing the superimposition of the radius artery waveform obtained through calcu- 
de^ce PrOCeSSm9 rad ' US 8rtery Wavetorm obtained throu 9h averaging processing in the pulsewave analysis 

Fig. 13 is a diagram for explaining the details of the normalizing processing which is applied to the radius artery 
waveform obtained as result of averaging processing in the pulsewave analysis device. 

Fig. 14 .s a diagram showing the correlation between the time f s of increasing pressure in the left cardiac ventricle 
and the time t p of one beat. 

Fig. 15 is a block diagram showing the structure of a pulsewave analysis device according to another embodiment 
of the present invention. 

H !f iS a . dia9ram showin 9 the conditions under which measurements employing a pulsewave detector 1 were 
conducted in this embodiment. 

Fig. 17 is a diagram showing the display of the superimposition of the measured and calculated waveforms of the 
radius artery, which is displayed by output device 6. 

i 1 f« S f ! qUint ViSW ° f lhe at,acnment of a sen sor 12 to band 13 of a wristwatch 1 1 , wherein the compositional 
elements 10 of the pulse wave analyser, excluding sensor 12, are housed in wristwatch 11. 

, n IS' 1 9 '. S 3 SqUint ViSW ° f 106 attachme nt of a sensor 22 to the base of a finger, wherein the compositional elements 
1 0 of the pulse wave analyser, excluding sensor 22. are housed in wristwatch 11. 

F.g. 20 is a structural diagram showing the attachment of sensor 32 and the compositional elements 10 of the pulse 
wave analyzer, excluding sensor 32, to the upper arm via a cuff. 

Rg. 21 is a block diagram showing the structure of a sphygmomanometer according to the second embodiment of 
the present invention. 

Fig. 22 is a diagram showing the relationship between the pressure waveform at the aorta (dashed line) and the 
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radius artery waveform (solid line) in a Type I pulsewave. 

Fig. 23 is a diagram showing the relationship between the pressure waveform at the aorta (dashed line) and the 
radius artery waveform (solid line) in a Type II pulsewave. 

Fig. 24 is a diagram showing the relationship between the pressure waveform at the aorta (dashed line) and the 
5 radius artery waveform (solid line) in a Type III pulsewave. 

Fig. 25 is a diagram showing the relationship between resistance R c in blood vessels at the center of the arterial 
system and distortion d. 

Fig. 26 is a diagram showing the relationship between resistance R p in blood vessels at the periphery of the arterial 
system and distortion d. 

10 Fig. 27 is a diagram showing the relationship between inertia L from blood flow and distortion d. 
Fig. 28 is a diagram showing the relationship between compliance C and distortion d. 
Fig. 29 is a diagram for explaining the contraction phase area method. 

Fig. 30 is a diagram showing the structure of a sphygmomanometer employing the conventional technology 
Fig. 31 is a flow chart of the program for obtaining capacitance C c . 
75 Fig. 32 is a flow chart of another program for obtaining capacitance C c . 

PREFERRED EMBODIMENTS OF THE PRESENT INVENTION 

< First Embodiment) 

20 

A first embodiment of the present invention will now be explained with reference to the accompanying figures. 

Fig. 1 is a block diagram showing the composition of the sphygmomanometer employed in this embodiment. In this 
embodiment, circulatory state parameters for the arterial system in the human body are evaluated based on information 
obtained from the body using non-invasive sensors. The specific details of these circulatory state parameters will be 
25 explained later. 

In Fig. 1, a pulsewave detector 1 detects the pulsewave at the radius artery via a pressure sensor S2 which is 
attached to the wrist of a test subject as shown in Fig. 2. Additionally, pulsewave detector 1 detects the blood pressure 
of the test subject via a cuff S1 attached to the upper arm of the subject as shown in Fig. 2. Pulsewave detector 1 cor- 
rects the measured radius artery pulsewave for blood pressure, and outputs the result as an analog electric signal. The 
30 analog signals are input to an A/D (analog/digital) converter 3, and are converted to digital signals at a given sampling 
period. 

Stroke-volume measurer 2 is connected to cuff S1 as shown in Fig. 2, and measures the stroke volume, i.e., the 
amount of blood which flows out from the heart per beat, via cuff S1. The results of this measurement are output in the 
form of a digital signal as stroke-volume data. A device which carries out measurements using a contraction phase area 
35 method may be employed for this type of stroke-volume measurer 2. 

Microcomputer 4 houses a waveform memory for storing the pulse waveforms which are taken up from the A/D 
converter 3, and houses a temporary recording memory which is employed as an operational region. Microcomputer 4 
carries out the various processing noted below in accordance with commands which are input from keyboard 5, which 
is an input device, and outputs the results obtained from processing to output device 6. The processing steps which fol- 
40 low will be explained in greater detail under the explanation of the present inventions operation which follows later. 

1. Pulsewave readout processing, in which time series digital signals of the radius artery waveform obtained via 
A/D converter 3 are taken up in the waveform memory (omitted from figures) which housed in microcomputer 4. 

2. Averaging processing, in which the radius artery waveform taken up in waveform memory is averaged at each 
45 beat, and a radius artery waveform (hereinafter, referred to as "averaged waveform") corresponding to one beat is 

obtained. 

3. Uptake processing, in which stroke-volume data is taken up in the temporary recording memory inside micro- 
computer 4. 

4. Parameter calculation processing, in which an equation expressing the radius artery waveform corresponding to 
so one beat is obtained, and each of the parameters in an electric model which corresponds to the arterial system is 

calculated based on this equation. 

5. First output processing, in which the obtained parameters are output to output device 6 as circulatory state 
parameters. 

6. Second output processing, in which the pulse waveform at the proximal portion of the aorta is determined from 
55 the obtained parameters, the systolic pressure value, diastolic pressure value, and the heart's workload at the prox- 
imal portion of the aorta are calculated, and these results are output to output device 6. 

Output device 6 will now be explained in detail with reference to Fig. 1 In this figure, the numeral 61 indicates a 
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measured Wood pressure display, which displays systolic and diastolic pressure values actually measured and the aver- 
age Wood pressure value. 62 is a display for displaying the estimated blood pressure value at the center of the arterial 
system. D,sp ay 62 displays average blood pressure E 0I . systolic pressure E m ; and diastolic pressure E 0 at the center 
of the arterial system which are obtained as a result of processing which will be described below. 63 is an alarm display 
which is composed of a plurality of LEDs which are disposed in a horizontal row. These LEDs light up in response to a 
difference between the systolic pressure £ m ' at the center of the arterial system and the systolic pressure value actually 
T, ^JIT y " e " the difference between tn * to"™ and the latter is less than ±10 mmHg, then a green LED 
stating NORMAL" is illuminated. Conversely, when the difference between the former and the latter exceeds +10 
mmHg. then a red LED stating "CAUTION" is illuminated. 

r-JU" a £ ar * met * r displa y When capacitance C c . electrical resistance R c , inductance L. capacitance C. electrical 
resistance R p , time duration / s of rising pressure in the left cardiac ventricle, time duration t D of a single beat stroke 
volume SV. and workload W s are supplied from microcomputer 4. parameter display 64 displays these parameters 
These parameters will be explained in detail below. 67 is a CRT display for displaying a variety of waveforms such as 
the waveform at the radius artery, the pressure waveform at the left cardiac ventricle, the pressure waveform of the 
aorta and the like. 65 .s a printer which prints out on paper the waveforms displayed on CRT display 67 and the various 
data displayed on measured blood pressure display 61. estimated central blood pressure display 62. warning display 
63, and parameter display 64. when print command button 66 is depressed. 

The significance of a warning display on warning display 63 will now be explained 

As explained for Figs. 22 through 24 previously, there are three types of differences which may be noted between 
the systolic pressure values of the estimated pressure waveform of the aorta and the radius artery waveform When the 
pulse waveform is a Type I variety (Fig. 22), then it is very likely that the test subject is a healthy individual. On the other 
hand, test subjects demonstrating Type II or Type III pulse waveforms are frequently in poor health 

For example, the Type II waveform (Fig. 23) is caused by an anomaly the state of blood flow, and is highly likely in 
pa bents suffering from a mammary tumor, liver or kidney ailments, respiratory ailments, stomach or intestinal ailments 
inf lammation, or the like. The Type III waveform, on the other hand, is caused by an increase in tension in the blood ves- 
sel wal s. and is very likely in patients having liver or gall ailments, dermatological ailments, high blood pressure, or pain 

Accordingly, the present embodiment provides that a warning display be carried out by means of illuminating a red 
LED when there is believed to be an anomaly in the difference between systolic pressure values as described above 

in the preceding example, a diagnosis was made based on the difference between the systolic pressure of the 
pressure waveform at the aorta and the systolic pressure value of the waveform of the radius artery. However it is also 
acceptable to employ the difference in minimum or average blood pressure values, instead of the difference in systolic 
pressure values. Moreover, it is of course acceptable to carry out diagnosis using the differences in maximum, minimum 
and average blood pressure values together. 

The present embodiment newly employs a "lumped five parameter model" as an electrical model for the arterial 
system. There are vanous components which determine the behavior of the circulatory system in the human body. 
From among these, the component of aortic compliance has been added to the four parameters of inertia due to blood 
at the center of the arterial system, resistance (viscous resistance) in blood vessels at the center of the arterial system 
due to blood viscosity, compliance (viscoelasticity) of blood vessels at the periphery of the arterial system, and resist- 
ance (viscous resistance) in blood vessels at the periphery of the arterial system, which were employed in the lumped 
four parameter model disclosed in Japanese Patent Application Hei 6-205747 (Title: Device for Analyzing Pulsewaves) 
o comprise this lumped five parameter model. An electric circuit has been employed to model these parameters Addi- 
tionally, we note here that compliance is a quantity expressing the degree of pliability of the blood vessels 

Fig. 3(a) shows a circuit diagram for a lumped four parameter model, while Fig. 3(b) shows a circuit diagram for a 
lumped five parameter model. The relationship between the parameters and the elements making up the lumped five 
parameter model is as follows. 



Capacitance C c : 
Electrical resistance R c 

Inductance L: 
Capacitance C : 
Electrical resistance R n 



to 



aortic compliance (cm 5 /dyn) 
blood vessel resistance due 
(dyn • s/cm 5 ) 

inertia of blood at center of arterial system (dyn 



blood viscosity at the center of the arterial system 

, . J/cm 5 ) 

compliance of blood vessels at periphery of arterial system (cm 5 /dyn) 
blood vessel resistance at periphery of arterial system due to blood viscosity (dyn • s/ cm 5 ) 

Currents /. i p , i c , and / s . which are flowing through each part of the electrical circuit, correspond to blood flow 
cm /s). Current / is the blood flow at the aorta and current is / s the blood flow pumped out from the left cardiac ventricle 
Input voltage e corresponds to the pressure in the left cardiac ventricle (dyn/cm 2 ). while voltage v, corresponds to the 
pressure (dyn/cm ) of the proximal portion of the aorta. Terminal voltage v p of capacitance C corresponds to the pres- 
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sure (dyn/cm 2 ) at the radius artery. Further, diode D shown in Fig. 3(b) corresponds to the aortic valve. Diode D is on 
(valve open) during a period corresponding to contraction, and off (valve closed) during a period corresponding to 
expansion. 

As will be explained below, these five parameters are not calculated all at once in the present embodiment. Rather, 
5 the lumped four parameter model disclosed in the reference cited above is employed to calculate all parameters with 
the exception of capacitance C c , after which capacitance C c is determined. A theoretical explanation of the behavior of 
the lumped four parameter model shown in Fig. 3(a) will now be made. 

The following differential equation is established for the lumped four parameter model shown in Fig. 3(a). 

10 di 

v x = R c i + L— +v p (1) 



Current / in the above equation may be expressed as: 
20 Thus, equation (1) may be rewritten as follows: 

25 



As is conventionally known, the general solution for a second order constant coefficient ordinary differential equa- 
tion may be obtained by summing the particular solution (steady-state solution) which satisfies equation (3) and the 
30 transient solution which satisfies the following differential equation. 

.-«>£ + (*o + i)5 + (i + £)i- (4) 

35 



The solution for differential equation (4) is obtained as follows. First, the damped oscillating waveform expressed 
by the following equation is assumed as the solution for differential equation (4). 

40 

v p = exp {st) (5) 
Substituting equation (5) into equation (4), equation (4) may be rewritten as follows. 

{«V+(lV7+£)'+(l+jj})}*-° ( 6 > 



50 

Solving equation (6) for s yields: 

- (*cC + jfc) ±J(R C C + ±) 2 - 4LC (l + fe) (7 ) 
s _ ___ 
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(*C + A) <4Lc(l + £) ,8, 
then the radicand of the radical sign in the second term becomes negative, and s is as follows. 



w 

s — 



(R C C + ^) ± jfiEc (1 + fe) - {R.C ^ 



75 



2LC 



Where, a is the damping factor and <o is the angular frequency. 

20 



a ~ 2LC (10) 
L + FLpR c C 



25 2LCRp 



30 



Further, when 

35 



40 



2LC 1 ; 



A, = LC 



a 2 = i±M5 (1S) 



50 



3 - —*r 

then equations (10) and (11) above may be expressed as follows: 



(14) 



55 
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" = Si (15) 



-ft* 



u = .,_-<*- (16) 



By confirming the value of s in this way, a solution may be obtained which satisfies differential equation (4). Based 
on the preceding view, then, equation (5) may be employed as an equation approximating the damped oscillating com- 
ponent which is included in the response waveform of a the lumped four parameter model. 

Next, the pressure waveform at the proximal portion of the aorta will be modeled. In general, the pressure waveform 
at the proximal portion of the aorta demonstrates a shape as indicated by the heavy line in Fig. 4. In this figure, t p is the 

20 time duration of a one beat in a waveform, while t s is the time period during which pressure increases in the left cardiac 
ventricle. In a lumped four parameter model, this pressure waveform is approximated by the triangular wave shown in 
Fig. 5. When the amplitude and duration of the approximated waveform are indicated as E ot E m , t p , and t p1 , then the 
aortic pressure v 1 at an optional time t may be expressed by the following equation. Here, E Q is the diastolicpressure 
(blood pressure during expansion phase), E m is the pulse pressure, (E 0 + E m ) is the systolic pressure (blood pressure 

25 during contraction phase), t p is the time duration of one beat, and t p1 is the time duration from the rise in aortic pressure 
until blood pressure falls to a minimum value. In the interval 0 < t < t p1 : 



0-£) 



v x = E 0 + E m 1 - — (17) 



In the interval t p1 < t < t p : 

35 v ^ = E o (18) 

The response waveform v p (i.e., radius artery waveform) when voltage v 1 which is expressed using equations (17) and 
(18) is input into the equivalent circuit shown in Fig. 3(a) is as follows. 
Intheinterval0<r< t p1 : 



v p = Enin + B (l - ^ + £> ml exp (-erf) sin (art 4-50 (19) 
45 In the interval t p1 < t < t p : 

^ = ^ + ^2* ex P {-« (f - fpf )} • sin {a> (r - t p1 ) + G 2 } (20) 

Here, E mjn is the diastolic pressure value in the radius artery waveform which is measured by pulsewave detector 
so 1 (refer to Fig. 1 1 explained below). 

The third term from the right in equation (19) and the second term from the right in equation (20) are the damped 
oscillating components in equation (5) explained above, a and co in these terms may be obtained from equations (15) 
and (16). 6, t b . D m1 , and D m2 are constants which are calculated in accordance with a procedure described below. 
Next, an examination will be made of the constants in equations (19) and (20), with the exception of a and co which 
55 were already confirmed. To begin with, the following equation is obtained when equations (17) and (1 9) are substituted 
into differential equation (3). 
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75 



20 



25 



30 



35 



40 



45 



w tb ( c Rp) 



|LC (* 2 - u, 2 ) D ml - Q D ml (r c C + -j^j + D m , (\ + j 



X exp (— at) sin (u;t + #i) 

+ 



ju>Z? ml ^i? c C + - 2LCauD ml ^ exp (-at) cos (ui + #i) 



The following conditions are necessary in order to establish the preceding equation (21). 

E 0 + E m = (^i + ^( Emin + B)=E 0 + A 3 B-f b A 2 (22) 



(21) 



E, 
t 



?m - B _i_ R '\ A * B , , 



(24) 



RcC + ~ = 2LCa (25) 
Hp 

so Equations (24) and (25) restrict a and co, however, these equations are satisfied by a and <d as obtained according 
to equations (15) and (16) above. 

When equations (18) and (20) are substituted into differential equation (3), the following equation is obtained. 



55 
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10 



15 



20 



25 



30 



35 



40 



45 



+ |LC (a 2 - a; 2 ) D m2 - a (fl c C + ^ £> m2 + (l + ^) D m2 j 
X exp {-a-(t - t P i)>sin{u (t - t,i) + #2} 

exp {-a{t - cos {w (f - i p i) + #2} 

(26) 

In addition to equations (24) and (25), it is also necessary to set up the following equation in order to establish 
equation (26). 



-H8 



E min = AsErnin (27) 



Next, each of the constants in equations (19) and (20) will be calculated based on conditional equations (22) 
through (25) and (27) for setting up differential equation (3). 

First, the following equation is obtained for E min from equation (27). 

E ■ = — ^ (28) 

^3 

An equation for B is obtained from equation (23): 

B = Mm (29) 

*pl^3 



When equation (29) is substituted into equation (22), and solved for t b , the following equation is obtained: 



= ^1^3 + ^2 (30) 



A 3 

50 Next, values are selected for the remaining constants D m1 , D m2% 0 7> and e 2 so that radius artery waveform v 

remains continuous over f=0, t p1 , t p , i.e. values are selected which satisfy the following conditions a through d. 

a * v P (tpi) in equation (19) and v p (t p1 ) in equation (20) are equivalent 
b. Vp(tp) in equation (20) and v p (0) in equation (19) are equivalent 
55 c. The differential coefficients at t-t p1 in equations (19) and (20) are equivalent 

d. The differential coefficient of equation (19) at f=0 and the differential coefficient of equation (20) at t=t are 
equivalent p 
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In other words, values are selected for D m1 and 6, so that: 



^ml 



(31) 



10 



15 

Where, 



0i = tan" 1 §i (32) 



(33) 



20 



D 12 = ( Voi -B- B min ) a + i + i£± (34) 



25 



and v o7 and i 01 are the initial values of v p and / c at /=0. 
Further, values are selected for D m2 and v 2 so thai: 



Dm2 = M±^ (35) 



35 



0 2 = tan" 1 ^ (36) 



40 

Where, 



D 2l =(^02 " E /77/J "© (37) 

45 



D 22 = (v 02 - E min ) -a + ^ (38) 

and v o2 and / o2 are the initial values of v p and / c at t=t p1 . 

In this manner, then, each of the constants in equations (19) and (20) are obtained. 

Next, the following equation for blood vessel resistance R c is obtained by means of an inverse operation from angu- 
lar frequency co in equation (16). 

55 
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5 



W 



20 



L - IFLpd LC (1 - u> 2 LC) 

Rc = CR, (39) 

Here, the condition for R c to be a real and positive number is: 

-*&-*L** HO, 



1 + (2u; J R p C) 2 o; 2 C 



15 in general, ft p and C are on the order of 10 3 [dyn • s/cm 5 ] and 10* 4 [cm 5 /dyn], respectively. a> may be viewed to be 
on the order of 10 (rad/s) or greater, since it is the angular frequency of the oscillation component which is superim- 
posed on the pulsewave. For this reason, the lower limit of equation (40) is viewed to be around 1/(a> 2 C). Accordingly, 
when L is approximated by the following equation (41) for the purposes of simplification 



R c becomes: 

25 

*c=^ (42) 

30 

From the relationship between equations (41) and (42), the damping constant a in equation (1 5) becomes: 

Q = rk (43) 

35 OKp 



Using the relationships between equations (41) through (43), the remaining parameters in the lumped four param- 
eter model may be expressed using a, a>, and L as follows. 

40 

R c = aL (44) 



45 



SO 



55 



^ = <*k (45) 



a 



C = 4y (46) 

Thus, it is clear that the parameters are confirmed by obtaining a, a> and L from the preceding equations (44) through 
(46). 
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a, co, B and t b may be obtained from the actual measured waveform at the radius artery and L can be calculated 
based on the stroke volume SV, as will be explained below The process for calculating L based on stroke volume SV 
will now be explained. 

First, the average value E 01 of the pressure wave at the proximal portion of the aorta is obtained from the following 
equation. a 

E 0l = 2 (47) 

tp 



The following equation may be established between R c , R p a, <a, and L. 

i?c + Rp = aL + ^ = (a 2 + J 1 ) ^ (48) 

The result of dividing the average current i.e.. average value E 01 , flowing through the lumped four parameter model 
t>y (H c +« p ) corresponds to the average value (SV/r p ) of blood flow flowing through the arteries due to the pulse 
Accordingly, the following equation may be established. 



SV 



{^+^)LT J) \ Eot ^-r~) < 49 > 



By solving the thus-obtained equation 49 for L, an equation for obtaining L from stroke volume SV may be obtained 
as follows. 



E 0 U + ^ 



L = a- op — ? — (50) 

■ ^ v « lue ^"-esponding to average current (1/f p ){£ 0 f E m IZ)} in equation (49) may be obtained by meas- 
uring the blood flow volume, and inductance L may be calculated based on this result. As a device for measuring blood 
flow volume, there are available devices employing the inductance method, the Doppler method and the like In the case 
of a device which measures blood rate using the Doppler method, there are available devices employing supersonic 
waves and laser. 

Next, a theoretical explanation will be made of the method for calculating the circulatory state parameters based on 
is a lumped f.ve parameter model. As mentioned before, the circulatory state parameters R c , R pi C, and L are determined 
using a lumped four parameter model. Accordingly, the value of capacitance C c is determined based on these param- 
eters. It is therefore necessary to obtain current ;. current / s , voltage v,, and voltage v p in Fig. 3(b). 

First the pressure waveform at the left cardiac ventricle is approximated with a sinusoidal wave such as shown in 
Fig. 4. In other words, by setting a> s =n/f s , the pressure waveform e of the left cardiac ventricle is expressed by the 
so following equation. 

e = E' sin co 



m s,n< °s r (51) 



Here E m ' is the systolic pressure and, if stated in terms of Fig. 5, corresponds to (E m +E 0 ). 

An explanation will now be made, taking the cases where time duration t corresponds to the contraction phase of 
f 'u" l*! 2 and t0 the expansion P hase fp+f,>. as shown in Fig. 4. Here, time f, and time u are the times at 

wh.ch the pressure waveform of the left cardiac ventricle waveform and the pressure waveform of the aorta intersect 
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(Contraction phase) 

In this case, the equation i/ r = e, as well as equations (1) and (2) for voltage v 1 and current /', respectively, are 
established. Accordingly, the following differential equation is set up from equations (1) through (3), (12) through (14) 
5 and (51). 

d v du 

Al ~d? +A2 ~dt + A * Vp = 3» 8inu; ** ( 52 ) 

w 

First, in the same manner as the lumped four parameter model, the steady-state solution v pst for this differential 
equation is obtained. For this purpose, the following equation is assumed for the steady-state solution v pst . 

15 

v p$t = E lC os<o s / + E 2 sinco s ? (53) 

By comparing the coefficients by substituting equation (53) for the term v p in equation (52), the following two equa- 
tions are obtained. 

20 

[A^JlAx) E x +uj 3 A 2 E 2 = 0 (54) 



25 



30 



35 



40 



45 



~ w,A 2 Ex + (A 3 - u*A x ) E 2 = E' m (55) 



Solving these equations, the following equations (56) and (57) are obtained. 

_ -u s A 2 E , rn 

El - (u s A 2 f + {A 3 -u*Ai) 2 (56) 



E 7 = {A 3 - JjAJ E* m 

(^sA 2 f + (A 3 -^A l ) 2 K } 



Next, the transient solution v ptr for differential equation (52) is obtained. Setting v p{r =exp(Xt) , this is substituted 
for v p in the following equation. 

50 
55 

As a result, the following equation is obtained. 

A,X 2 + A 2 X + A 3 = 0 (59) 
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Solving this equation for X, the following equation is obtained. 



A _ -A 2 ± y /A 2 2 -AA l A z = / A 2 \ + / /j 2 V 
" 2,4, ~V 2Aj~ \\2aJ Ai 



( a 2 y _ A3 



(60) 



70 



75 



20 



Setting {4^/(24 7 )} 2 <(A 3 /A 1 ) (i.e., oscillation mode), the following equation is obtained. 



A = 



A 2 



In this case. 



A * A - 



2Aj 



(61) 



(62) 



25 



U/i 



(63) 



30 Here, transient solution v ptr is set as in the following equation. 

v ptr = ( a 1 cosco 1 f + ya 2 sina> t t) exp (-p 1 f) 



(64) 



As a result voltage v p is expressed as the sum of the steady-state solution and the transient solution, and therefore 
35 may be obtained from the following equation using equations (53) and (64). 



40 



45 



50 



v p = (E, cosco s f + E 2 sina> s f) + ( ttl coscd^ + ja 2 sino> 1 1) exp (-p 1 1) 
The following equation may be obtained for current / by substituting equation (65) into equation (2). 



-=r I sinuj.t 



- + u; 3 Ci? 2 ) cos a,,* + (-u^Z?, + ^ 



(66) 



(65) 



55 



Next, the v 02 and , 0 are assumed as the values of v p and /. respectively, when t=t , as in the fallowings equations. 

'0 = J 0 + (<*,,;, +ya 2 J 2 )exp(-p 1 / 1 ) (67) 
v 02 =P 0 + (a,P 1 +/a 2 P 2 )exp(-p lfl ) {68) 
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As a result, equations (65) through (68) may be used to establish the following equations. 



Jo = 



- (j^ + u s CE^j cosoj a t x + ^-u,CEi 4- y£j sin u: 3 ti 



( I -Pi 



J2 = U>iC COSUJit l -h 



(HP) 



(69) 



(70) 



(71) 



20 



25 



30 



P Q = C0Stt> s ^ + E 2 sinoo^ (72) 

P n = cosco^ (73) 

P 2 =sinc0 1 f 1 (74) 
Next, solving equations (67) and (68) for a 1 and a 2 , the following equations (75) and (76) may be obtained. 

_ f (^02 ~ Pp) J 2 ~ (jo ~ Jo) P 2 



ai = 



exp (/?iia) 



(75) 



35 



a, - / -fote-^i + fa- Jo) Pi ) ,„ . . 



(76) 



40 



45 



50 



55 



Next, the relationship expressed by the following equation may be established from equations (70) through (74). 

(77) 



J 2 P 1 ' J 1 P 2 =CD 1 C 



Accordingly, when equations (75) and (76) are substituted into equation (64), and equation (77) is employed, then 



the following equation is obtained as transient solution v t 



Vptr = 



ptr 



(^02 - Po) COSCJ! (£ — t x ) — 

{(1 - PiCRp) {vq2 - Pp) - Rp (ip - Jp)} {sinu^ (t - fr)} 
x exp exp (-/?!*) 



(78) 



Where, when 
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b Mr = V Q2 ' H Q 
V = t - f , 



B 2tT = JJ^l£i£3z11^^ 

^"iCR, (81) 

then the following equation (82) is obtained. 

v P tr = ( a i fr cosa> ,\ ' + B 2 tr sin a> , V ) exp (-p 1 f ) 
Accordingly, equation (65) becomes as follows. 

v p = (E 1 cos© s f + E 2 sin<o s 0 + (S 1(r cosa> ,r + B 2f ,sin co 1 f) exp (-p t r) 
Next, the tennis in the above equation (66) are defined as follows: 

_ E 1 

D ist = — 4- u s CE 2 (84) 



27» 

D 2st = -uj s CE 1 + ~ ( 85 ) 

Dur = ^^^ Bitr + UiCB* (86) 



»r (87) 



As a result, the following equation (88) is obtained for current /: 

' = ( D i5f cos< » s f + D 2s ^inw s t) + (D n ^cose^ f + D 2u s\n<i>,r ) exp (-0 ,f ) 
Current / s is obtained as the following equation. 
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(Expansion phase) 

During expansion, diode D turns OFF, ending the impression of the left cardiac ventricle pressure e on the cathode 
side of the diode D circuit. As a result, the current flowing through capacitance C c has the same absolute value and the 
opposite polarity of current /. Accordingly, voltage ^ may be expressed by the above equation (1), while currents / and 
i c may be expressed by the following equations, respectively. 

i = -C C £ ,90) 



Accordingly, voltage v p becomes 

*-H,(«-<0-*(ftil + C*k) (92) 



Further, since /-/c= i p -v p lR pt the following equation results. 

, = ^ +C ^ (93) 



35 Substituting equation (93) into equation (1), and differentiating both sides of the obtained equation with respect to 
time t, the following equation is obtained. 



40 



dt ' dt 3 



45 



SO 



55 



Next, the following equation may be drawn from equations (90) and (93). 

dvi_ % + C% 
dt c c 



Further, the following equation may be obtained from equations (94) and (95). 



(95) 



LC n*+{ r, )^+{ — cX )-t + {c^h =0 (96) 
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Accordingly, this equation may be restated as: 
d?v p , d?v p A . dv v 

■5r + A-^ + ^+^tv = o (97) 



where: 

w 



15 



20 



25 



30 



A , _ L + CRcR, 
1 LCR^~ (98) 



A , CcR; + cjt, CIL, 
2 LC^r% " (99) 



4 ~ Zqc^ ( 10 °) 



Next, the following equation is obtained when v p =exp{\t) and substituted into equation (97). 

(X 3 + A\X Z + A' 2 X+ A' 3 )exp{Xt)=0 ( 101 ) 
35 Next, the following definitions are provided. 

._ A' t 2 A' 2 

P ~ ~9~~T (102) 

40 

Ai 3 A 1 A 1 A 1 
9 27 + — ~ T ( 103 ) 

45 



r~2. — 3 1/3 



(104) 
(105) 



55 



20 
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cl = - {u + v) + ^. 



(106) 



02 = 



u + v A\ 

— + i 



(107) 



10 



U>2 = 



(u-v) — 



(108) 



15 



1 = (109) 

20 x 2 =-p 2 +yo> 2 (no) 

A, 3 = -p 2 -ya> 2 (111) 

Note that the oscillation mode is indicated when (cf-p 3 ) > 0. 
25 Next, voltage v p is assumed according to the following equation. 

v p = b, exp(-a-f) + b 2 exp{(-p 2 + ja> 2 )f} + b 3 exp {(-p 2 - j<x> 2 )t} (112) 

As a result, current / may be expressed by the following equation after substituting equation (112) into equation 

30 (93). 

/ = g 0 6 1 exp(-a'f) + (g, +jg 2 ) b 2 exp {(-p 2 + jo> 2 )t } + (g r jg 2 )b 3 exp {(-p 2 - j(o 2 )t} (113) 

where, 

35 

9o 

40 

9\ 

45 

92 

50 

Accordingly, voltage v, becomes as follows from equation (113). 



1 - a'C/Z, 
1 -/? 2 Cfip 



^CRp 



(114) 



(115) 
(116) 
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10 

where, 



15 



Vi = -k J* * 

= /o6i exp(-a't) + (/ 2 + jf 2 ) ^ exp {(-/3 2 + ju 2 ) t} 

+ (/i - jh) h exp {(-/? 2 - ju> 2 ) t} p 17 j 



20 



25 



71 " C^ + u,§)C7 c ("9) 



f„ _ ^2^1 + 0292 

<#+«!)C e (120) 



Next, for the purpose of simplifying calculations, time t 2 shown in Fig. 4 will be set to t=0. When voltage v, voltage 
so v ? and current ; at f=0 are set respectively to v 01 , v 02 and i 0 , then the following expressions are obtained when the term 
t in equations (1 1 7), (1 12) and (1 1 3) is set to f=0. 

v oi = ;C o 6 i +d +jf z )b 2 +(f^ -jf 2 )b 3 (121) 

v 02 = t> 1 + 6 2 + o 3 ( 12 2) 

'o = 0o^i + y0 2 )f> 2 + (0 1 -yg 2 )£> 3 (123) 
Rewriting the second and third terms in equation (1 12), the following equation is obtained for voltage v„ 

40 " 

v p = exp(-a'f) + {(6 2 + fo 3 ) cos<o 2 f + j {b 2 - 6 3 ) sina> 2 f} exp (-p 2 f) (124) 
= B 0 exp (-a t) + (B 1 cos© 2 1 + B 2 sin<» 2 1) exp(- p 2 f) 



45 Where. 



^452 — £3/2 

50 



55 
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5 
10 

and 

k ^ = v o^ ' f o v 02 (128) 

15 

k 2 = 'o " 00^02 (129) 
^3 = ^1*^0 030) 
*° * 4 -'i-'o (131) 

Next, the following equation is obtained for current / by substituting v p of equation (124) into equation 2. 

/= D 0 exp(-a'f) + (D, cos<x> 2 t+ D 2 s\no> 2 t) exp (-p 2 f) (132) 

25 

Where, 



30 



35 



40 



45 



D ° - (HP)* 



(133) 



£» 2 = -u) 2 CB x + ( lz^C^ \ 



£0 



Accordingly, from equation (132), voltage ^ becomes: 

= H 0 exp(-a't) + (Hi cosuj 2 t + H 2 sinu/ 2 £) exp (-/? 2 t) (136) 



55 



Where, 
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H Q = 

5 
10 

H, 



Do 

a'Cc (137) 
{0l + ul)C c (138) 



= ~uj 2 D x +(3 2 D 2 



75 



20 



m Jh 1° P reCedin , 9 e f P' a nation. time t 2 was set to r=0. Therefore, a substitution of t -» (f- r2 ) is carried out in order to 
match the time scale. Voltage v,. voltage v p and current / can be obtained as follows from equations (136), (124) and 
respectively. x ' 



H = J^x eX ^ { ; a 7V/S )} {W i cos «>2 (' - r 2 ) + « 2 sinco 2 (f - f 2 )} exp {-p 2 (f - f 2 )} (140) 
M„ exp {-a (r - f 2 )} [H m sin [ &2 (f - f 2 ) + + 2 ,}] exp {-p 2 (t - t 2 )} 

25 V PJ B o {-«* (' - ' 8 H + {S, cos© 2 it - f 2 ) + e 2 sina) 2 (f - f 2 )J exp {-p 2 (f - f 9 )l (141) 

= e 0 exp { -a' (t - t 2 )} + [B m sin { <o 2 (t - t 2 ) + + 2 2 }]exp {-p 2 (f - t 2 )) 2 1 ' 

/ = D Q exp {-«' (f - t 2 )) + (D, cose, (r - f 2 ) + D 2 sin<D 2 (f - t 2 )} exp {-p 2 (f - t 2 )} (142) 

= D 0 exp {-« {t - t 2 )} + [D m sin {<* 2 (t - t 2 ) + $ 2 3 }]exp {-p 2 (f - f 2 ) } 2 " 1 ' 



30 



35 



40 



45 



Where, 



(143) 



<*>2i = tan- 1 ^- 

# 2 (144) 



(145) 



<f> 22 = tan- 1 ^! 

£ 2 (146) 



<t>23 = tan" 1 ^ 

D 2 (148) 



(147) 



Here, current / s during the expansion phase is 0. 

Next, a theoretical value for stroke volume SV will be obtained. Stroke volume SV is given by the area of current / 
during the contract.on phase, and thus may be obtained by integrating equation (89) for current / s over the range from 
tf to i2- Namely: 
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SV = f 7 i s dt 




5 



■cB' m + D ut \ , . 
(si 





(cOSU s t 2 — COS LJ S ti) 



10 



+ 



75 



+ 



(149) 



20 



An explanation of the operation of the pulsewave analysis device according to the present embodiment will now be 
made. Figs. 6 through 10 are flow charts showing the operation of this pulsewave analysis device. Fig. 1 1 is a waveform 
diagram of the average waveform which is obtained as a result of the averaging processing described above. Fig. 12 is 
a waveform diagram contrasting a radius artery waveform obtained by parameter calculation processing, to be 
25 described below, and the average waveform which is obtained as a result of averaging. An explanation of the pulsewave 
analysis device's operation now follows, with reference given to these figures. 

<3) Pulsewave readout processing 

30 In order to evaluate the circulatory state parameters in a test subject, an diagnostician attaches the cuff Si and 
pressure sensor S2 to the test subject as shown in Fig. 2, and inputs a start-measurement command via keyboard 5. 
In response to this command, microcomputer 4 sends an indication to pulsewave detector 1 to measure the pulsewave. 
As a result, pulsewave detector 1 detects the radius artery pulsewave, and A/D converter 3 outputs a time-series digital 
signals showing this radius artery pulsewave. Microcomputer 4 takes up the digital signals in the waveform memory 

35 housed therein over a fixed interval of time (approximately 1 minute). In this manner, the radius artery waveform for a 
plurality of beats is taken up in waveform memory. 

@ Averaging processing 

40 Next, at each beat, microcomputer 4 superimposes the radius artery waveforms over a plurality of beats, and 
obtains the average waveform per beat over the aforementioned fixed interval of time. Then, this average waveform is 
stored in the internal memory as a representative waveform of the radius artery waveforms (Step S1). An example of a 
thus-produced representative waveform W1 of the average waveform is shown in Fig. 11. 

45 ® Uptake of stroke-volume data processing 

Next, microcomputer 4 sends an indicator to measure stroke volume to stroke-volume measurer 2. As a result, 
stroke-volume measurer 2 measures the stroke volume in the test subject, with the measured result taken up by tem- 
porary memory in microcomputer 4 (Step S2). 

50 

@ Parameter calculation processing 

First, the four circulatory state parameter excluding capacitance C c are determined based on the lumped four 
parameter model. 

55 Processing in microcomputer 4 proceeds to Step S3, executing the parameter calculation processing routine 
shown in Figs. 7 and 8. At this time, the routine shown in Fig. 9 for calculating a and co is carried out (Steps S109, S1 1 7). 
Accompanying the execution of the routine for calculating a and a>, the routine for calculating a> shown in Fig. 10 is car- 
ried out (Step S203.) 
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An explanation w,ll now be made of the details of the processing carried out in the routines mentioned above First 
nfj^T 9 " * m,crocom P uter 4 determines for the average waveform of the radius artery the time t , ' and Wood 
pressure value y, corresponding to a first point Pi at which blood pressure reaches a maximum value; the time fe'and 
btood Pressure value y 2 corresponding to a second point at which blood pressure falls subsequent to the first point- the 
time fe and blood pressure value y 3 corresponding to a third point P3 which is the second peak point and the time t 
per beat and diastolic pressure value E min (corresponding to the first term in equations (3) and (4) above) (Step S10lf 

' nJH li °T » Where ?i S diffi ° Ult t0 discriminate between second point P2 and third point P3 when the pulsewave is 
gentie. then the times at the second point and the third point are assumed to be f~' = 2?/ and /-'-a/ respectively. 

Next, in order to simplify the processing, normalization of blood pressure values y, to y, is carried out usinq Wood 
pressure value y 0 at point A shown in Fig. 13 (Steps S102. S103). and the value at point B is initially set to (y^l-O 1 
Next, optimal values for B. t b< a and © are determined according to the following process. ' 

a) B is varied within the range l(y</2)~y 0 ] and t b is varied within the range l(tp/2)~t p ], at +0.1 intervals in both 
cases. The a and o for which | v^-y, 1. 1 Vf /t 2 7-y 2 1 and | v p (t 3 )-y 3 1 are minimized are determined for each B and 

b) From among the B, tb . a and «, values obtained in a) above, the B. t b a and «, for which | v p (t,-)-y ,|, I v D fM-y P 
| and | v p (t 3 ')-y 3 | are minimized are determined. p 1i ' pl 21 Ys 

c) Using the B f and t b obtained in b) above as standards, the processing in a) and b) is carried out again for B within 
•n the range of S+0.05 and for within the range of fct0.05. 

d) When carrying out the processing in a) through c) above, a is varied by 0.1 increments within the range from 3 

r<J « T/^T Z ^ Z 6 Ca,CU,ated for each «■ Next ' for *«* a, dichotomy is used to obtain the <* at which 
dv (t 2 )/dt=0 (see the flowchart in Fig. 10). When calculating the value of v p in the above processing the initial 

are U deter'n!inS Uat '° n '* ** tC> ^ * ^ * preceding Processing, the final values of B. t b a and a> 

e) fp,. E m and E 0 are calculated based on equations (28) through (30) and (44) through (46) (Steps S123 S124) 

f) The value of L ,s calculated based on the measured stroke volume, using equation (50) (Step S125) and the 
remaining parameters R c , R p . and C are obtained from equations (44) through (46) (Step S126). 

thf' baS6d 00 ' Umped f,Ve P arameter model, the final circulatory parameter of capacitance C c is determined 

thJl^h IT 13036 ' Hf*? 8 T 56 considered wher ein capacitance C c is determined so that the calculated value and 

catcSS S? 6 / T° ke V °!T e SV ^ 6quivalent and Wherein capacitance C c is determined so that the 

mi?hS^Th , mea , $Ured Va ' Ue 01 thS diaStoMc pressure of the P"'sewave are equivalent. Each of these 

methods will be explained separately. 

®-1 . Method to determine capacitance C c (aortic compliance) so that the calculated and measured values of stroke 
volume are equivalent 

m^nr!H Pl f. nati °r r"J irSt !* m o de ° f 3 Sp6Cif iC meth0d ,or determini r>g capacitance C c so that the calculated and 
measured value of stroke volume SV are equivalent. 

whvT' 81 ' the 1 value ° f caprice c c is estimated as in the following equation, based on the value of capacitance C 
wh,ch was calculated using the lumped four parameter model. Moreover, values obtained using the lumped four param- 
eter model are employed for the values of the other circulatory state parameters, R^ R p , C and L. 

c c = 10 * C (150) 

M « N .!f"+ tr01 ? ^ IUme SV iS calcu,ated according to equation (149). using these circulatory state parameters. In this 
case the time f s of nsing pressure in the left cardiac ventricle is estimated according to the following equation from the 
time t p of one beat which was obtained from the lumped four parameter model. 

t s = (1.52-1.079f p )j p (151) 

This relational equation is an experimental equation obtained as a result of measuring the contraction time of the 
left cardiac ventricle using a heart echo. As shown in Fig. 14, 

-0.882 is obtained as the correlation coefficient. Additionally, the value obtained using the lumped four parameter model 
I J^rl ° r t systolic pressure E -' < see *V*ions (22) and (28)). In addition, time f, and t 2 are obtained from the 
relationship: 'eft cardiac ventricle pressure=aortic pressure. Moreover, since v 02 and i 0 are the values of v D and / at 
r=/,.f, is substituted forfm equations (83) and (88), to obtain v 02 and i 0 . P 

The value of capacitance C c is determined so that the calculated value for stroke volume SV obtained as above is 
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equivalent to the measured value taken up from stroke-volume measurer 2. In other words, the value of capacitance C c 
is varied within in a fixed range starting from an initial value obtained from equation (150). Then, the measured value of 
stroke volume and the value which was calculated from each capacitance C c values is compared, and a check is made 
to see if the integer part of the measured and calculated values are equivalent. If the integer parts are equivalent, then 
5 the measured and calculated values of capacitance are deemed equal, capacitance C c is determined, and parameter 
calculation processing is terminated. 

On the other hand, if no coincidence is observed in the measured and calculated values of stroke volume merely 
by adjusting the value of capacitance C c ,then from among the adjusted capacitance C c values, the capacitance C c 
value at which the difference between the measured and calculated values for stroke volume is minimized is selected 
io as the final value for capacitance C c . Next, the value of systolic pressure E m ' is varied by 1 mmHg within a range of ±3 
mmHg, and a check is made for the presence of a coincidence between the calculated and measured values for stroke 
volume in the same manner as described above. If a systolic pressure value E m ' is present at which coincidence is 
observed, then that value is set as the final systolic pressure value E m ' , and parameter calculation processing is termi- 
nated. 

75 If no coincidence is observed between the measured and calculated values for stroke volume even after adjusting 

the value of systolic pressure E m ' , then the value of resistance R p is adjusted. Then, from among the adjusted systolic 
pressure values E m ' , the value at which the difference in the measured and calculated values for stroke volume is min- 
imized is set as the final systolic pressure value E m ' . Next, resistance R p is increased or decreased by intervals of 
1 0fdyn • s/cm 5 ], for example, with the R p value at which the difference between the measured and calculated values for 

20 stroke volume is minimized selected as the final resistance R p value. 

An example of a flow chart for realizing the above-described process is shown in Fig. 31 . Note that with respect to 
the parameters which are varied within a fixed range in the program, a subscript V has been added to indicate the orig- 
inal parameters. 

25 <g)-2. Method for determining capacitance C c so that the diastolic pressures of the calculated pulsewave and the meas- 
ured pulsewave are equivalent 

Next, an explanation will be made of a method for determining capacitance C c so that the diastolic pressures of the 
calculated pulsewave and the actual measured pulsewave are equivalent 
30 In this method, a contraction period QT is first obtained in advance from the test subject s electrocardiogram. Next, 
with respect to this contraction period QT, the time t s v during which pressure in the left cardiac ventricle is rising is var- 
ied by intervals of 0.01 sees within the range from [QT+0.1(sec)] to [QT+0.2(sec)]. At the same time, systolic pressure 
E m V is varied by intervals of 1 mmHg over the range from [E 0 +E m -20 (mmHg) ] to [E 0 +E m +20 (mmHg) ). 

In other words, the total number of combinations of each value of time t s v of rising pressure in the left cardiac ven- 
35 tricle and systolic pressure E m v' is 451 . A capacitance C c is calculated for these combinations so that the diastolic pres- 
sure values of the calculated pulsewave and the measured pulsewave are equivalent. 

Next, the sampling value for the calculated pulsewave in each combination is defined as P^t) and the sampling 
value for the actual measured pulsewave in each combination is defined as P 2 (t), and the average deviation e of the 
waveform in each combination is obtained from the following equation. The capacitance C c (aortic compliance) at which 
40 this waveform average deviation e is minimized is employed. An example of a flowchart for realizing the above- 
described process is shown in Fig. 32. 

£ =|!^)^ (152) 



In this manner, then, the circulatory state parameters at which the measured and calculated values of stroke vol- 
50 ume are equivalent are all determined. 

© First output processing 

Once the above-described process for calculating the parameters is completed, microcomputer 4 outputs the cir- 
55 culatory state parameters L. C. C Cj R c> and R p to sequential output device 6 (Step S4). 

The values of circulatory state parameters calculated from the radius artery waveform in the case where the test 
subject is a 32 year old male are shown below. 
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Capacitance C c 


= 0.001213 [cnrVdyn] 


Electric resistance R c 


= 98.768 [dyn-s/cm 5 ] 


Inductance L 


= 1 5.930 [dyn-s^cm 5 ] 


capacitance o 


= 0.0001241 [cm 5 /dyn] 


Electric resistance R p 


= 1 300.058 [dyn-s/cm 5 ] 


Time t s of rising pressure in left cardiac ventricle 


= 0.496 [s] 


Time t p of one beat 


= 0.896 [s] 


Stroke volume SV 


= 83.6 [cc/beat] 


Systolic pressure E m ' ^ j 


= 117.44 [mmHg] 



.,.i a ^ ShOW ? R ? J 2, thera iS 900d coincidence between the measured waveform of the radius artery and the cal- 
culated waveform of the radius artery obtained from the calculated parameters. 

© Second output processing 

L r n r te Z S !' * he P re K f ure , w f eform °* ^e aorta is obtained based on the values of the circulatory state parameters 
°' ° c ' * c ' and £»■ Nam e'V. by employing equation (51) during the contraction phase and equation (140) during the 

^rr^ir?^ WaV f rm ° f V ° na9e V1 iS calculated for one b *at only (ie.. time 0~time f p or time ,,~time ( f?+ U). 
Then, the calculated waveform .s output to output device 6. with the pressure waveform of the aorta displayed Next 
l^h ° ° btained waveform of the proximal portion of the aorta at time t, is calculated from these equations' 
and th,s calculated value is defined as the diastolic pressure value £ 0 . The thus obtained systolic pressure value E ' 
.s relayed to output device 6 together with the diastolic pressure value, and these values are displayed on output device 

« J^Z^ 9 * by means of *e Present embodiment, it is possible to display each of the circulatory state parameters, 

Itl 1 ^ ? S 1 T 6 ! 8 Ur !: d ' aSt0liC preSSUre ** * e pressure ^veform of the aorta at the center of the arterial 
system, for the test subject or diagnostician. 

f urthe ^ s 'nce equation (51 ) indicates the pressure waveform at the left cardiac ventricle, the pressure waveform of 
the left cardiac ventricle may be displayed on output device 6 in place of the pressure waveform of the aorta, as the 
pressure waveform at the center of the arterial system. 

(Embodiment 2 > 

In the preceding first embodiment, the values of the circulatory state parameters were calculated from the radius 
artery waveform and the stroke volume. However, as explained above, in order to carry out a detection of stoke volume 
l!!! CeSSary toattach ^ S1 to 106 test su °iect, subjecting the subject to some inconvenience. Accordingly this 
embodiment provdes for an estimation of blood pressure values at the center of the arterial system by focusing atten- 
-on on a change ,n aortic pressure based on the shape of the radius artery waveform, and using distortion to represent 
the shape of the waveform. In other words, in this embodiment, the circulatory state parameters are derived based on 
distortion d which is obtained from the radius artery waveform. 

First, in the same manner as in the first embodiment, microcomputer 4 carries out ® pulsewave readout and © 
averaging, to obtam the average waveform per beat for the radius artery waveform. Next, a Fourier analysis of the 
pulsewave is performed by carrying out a fast Fourier transform (FFT) on the average waveform. Next, the fundamental 
wave s amplitude A, , the second harmonic wave's amplitude A 2 , the third harmonic wave's amplitude A 3 to the nth 
harmonic wave's amplitude A n are obtained from the frequency spectrum obtained as a result of this analysis. Here the 
value of n (which ,s a natural number) is optimally determined after taking into consideration the size of the amplitude 
o Mhe harmonic waves. Based on these amplitude values, then, distortion d defined by the following equation is calcu- 
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distortion i= M±3±l^Wl (1M) 



A 



1 



Next, the circulatory state parameters are estimated from the obtained distortion d. This estimation is carried out 
based on the understanding that there is a correlation in the degree of correspondence between the distortion of the 

10 radius artery waveform and each of the values of the circulatory state parameters. Namely, distortion d and circulatory 
state parameters are measured in advance for a number of test subjects, and a relational equation between distortion 
and the circulatory state parameters is derived. Examples of correlations obtained as a result of measuring distortion 
and the circulatory state parameters R c , R p , L and C are shown in Figs. 25 through 28. Aortic compliance C c is not 
shown in these figures, however a correlation coefficient and a relational equation therefor may be obtained in the same 

is manner as for the other four parameters. 

Based on distortion d calculated from the above equation (153) and the relational equations shown in each of Figs. 
25 through 28, the circulatory state parameters of R c , R pt L, C, and C c are calculated. Next, output processing in steps 
© and <g> of the first embodiment are carried out in the same manner, to obtain the waveform of one beat of the pres- 
sure waveform of the aorta from the calculated circulatory state parameters. At the same time, diastolic pressure value 

20 E 0 and systolic pressure value E m at the proximal portion of the aorta are calculated and displayed on output device 6. 

(Embodiment 3 > 

In this embodiment, in addition to systolic and diastolic pressure values at the proximal portion of the aorta, the 
25 workload on the heart (hereinafter, referred to as "cardiac workload") is calculated from the blood pressure waveform 
at the proximal portion of the aorta which is obtained as described above, and is displayed. 

Cardiac workload, one indicator showing the load on the heart, is defined as the product of stroke volume and aor- 
tic pressure, and is the result of converting stroke volume per minute to workload. 

Here, stroke volume is defined as the volume of blood flow sent out from the heart at each pulse, and corresponds 
30 to the area of the waveform of blood flow from the heart. Stroke volume is correlated with the area of contraction phase 
of the pressure waveform of the aorta, and may be obtained by applying the contraction phase area method on the 
pressure waveform of the aorta. 

In other words, the area S under the portion of the pulse waveform corresponding to a contraction phase in the 
heart is calculated. Explained in terms of the pulse waveform in Fig. 29, this is the area of the region from where the 
35 pulsewave is rising to where it notches, as indicated by the hatching. Next, stroke volume SV is calculated using the 
following equation, where K indicates a specific constant. 

Stroke volume [ml] = area S [mmHg • s] x K 

40 Cardiac stroke volume is defined as the volume of blood flow sent out from the heart per minute. Accordingly, car- 
diac stroke volume can be obtained by converting stroke volume to minutes. In other words, cardiac stroke volume may 
be obtained by multiplying stroke volume by heart rate. 

In this embodiment, in the output processing in © of the first and second embodiments, microcomputer 4 calcu- 
lates the cardiac workload based on the pressure waveform calculated for the left cardiac ventricle, and displays this 

45 result on output device 6. The other processing is equivalent to the preceding embodiments, and an explanation thereof 
will thus be omitted. 

Microcomputer 4 calculates cardiac workload W s according to the process shown below. 

First, w s is defined as the product of e • i & and may be calculated as follows from equations (51), (88), and (89). 

50 . 2 

w s = e-/ s =a> s C c E' m sinot> s fcosa> s f + E' m sin<o s f (D 1sf cosa> s r + D 2sf sina> s f) + E' m sina> s f 

(D Ur cosa> t V + D 2 fr sincD 1 V ) exp (-p 1 f) 

55 Setting the first, second and third terms in equation (154) to w 1t w 2 and w 3 respectively, these terms may be 
rewritten as follows. 
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oj C E* ^ 

Wl = -^-f^sin^t (155) 



10 



E' 

U2 - ^{A 5i sin2c^-Z^ 4 ( cos 2a/ s t- 1)} 



E f 



(156) 



~A>tr {cos (u, a t + cj^') - cos (u s t - tJit')}] exp (-pit 9 ) (157) 

20 Next, using equation (80), the following expressions are defined: 

co s r + 00^' = co s f + (f - fl ) = (co s co 1 t 1 = © a f - ci> (158) 

25 <» s t-G> 1 f'=(«> s -co ,) f + to ^ =co b r + <D (159) 
Namely, 



30 



35 



40 



45 



50 



co a =co s + co ! (160) 
co b = « 5 - co 1 (161) 

<t> = e > 1*1 (162) 
In this case, then, equation (157) becomes as follows: 

™3 = -^Piir{sin(a; a t-*)+sin(a; 6 t + $)} 

~D 2tr {cos (u; a t - $) - cos [u b t + <£)}] exp [-0it') (163) 

/i«wf;^: ^JS? ^ 3 respective, y defined as fo,low s. and the following equation is derived from equations 
( \oo), (157) and (163). 
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/C E' 2 
w l dt = — ^-cos 2u s l 



CcE' 2 

c m 



(l - 2cos 2 u; s t) 



(164) 



70 



75 



W 2 = I W2 dt 

= Wist + 2D 2at u a t) 



2m 5 t 
2oj s 



-0} 



20 



- 2 cosa;^ cos^t + Z) 2 5t sina; s t)} 



(165) 



25 



30 



35 



40 



45 



/E 9 
w 3 dt = 
2 



(-U a D Ur +/3iD 2t r) cos {u a t - <i>) 
— (P\D ltT + oj a D 2t r) sin (a7 0 * - <£) 



(a; t £>i tr + piDztr) cos (a; fc < + $) 

+ {-/3iD Ur -+- uj h D 2t r) sin (u/ 6 i + $) 



/3f + w? 



► exp (-AO 



(166) 



50 Since workload W s is obtained by converting the sum of W u W 2 , and W 3 to a "per minute value," it may be 

expressed finally as the following equation. 
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The significance of displaying cardiac workload, in addition to systolic and diastolic pressure values at the proximal 
portion of the aorta, as explained above is as follows. 

Namely, in the conventional sphygmomanometer, blood pressure is measured at the periphery of the arterial sys- 
tem such as at the upper arm. radius artery or the like, and is therefore an indirect measurement of the load on the 
heart. However, changes in the load on the heart may not always be reflected in peripheral blood pressure values 
Accordingly, using the periphery of the arterial system to interpret load on the heart cannot be viewed as a method of 
absolute accuracy. 

Accordingly, in the present invention, emphasis is placed on the importance of using blood pressure waveforms at 
the center of the arterial system to view cardiac load, with the blood pressure waveform of the proximal portion of the 
aorta (i.e.. the center of the arterial system) estimated from pulse waveforms measured at the periphery. If the systolic 
and diastolic pressure values at the proximal portion of the aorta are then calculated from the estimated pressure wave- 
form at the aorta, then these blood pressure values become indicators which directly show the load on the heart 

By obta.n.ng cardiac workload as described above based on the pressure waveform at the aorta it is possible to 
provide thesystohc and diastolic pressure values at the proximal portion of the aorta as other useful indicators of car- 
diac load. Therefore, an explanation of the usefulness of cardiac load will now be explained using examples 

Fn-st the case will be considered in which a hypertensive agent is administered to a patient to treat high blood pres- 
sure. Ordinarily, if the drug is effective, a change will be noted in the systolic and diastolic pressure values which are 
measured at the radius artery. However, even if no change is observed in the systolic and diastolic pressure values it 
is still poss-ble that the agent is having an effect in easing the load on the heart. This is because it is not absolutely 
essential that a hypertensive agent reduce the blood pressure at the radius artery, but rather, the agent's function is suf- 
ficient if it reduces the load on the heart somewhere along the arterial system. Thus, even if no marked change is 
observed I in the blood pressure values at the periphery of the arterial system, such as at the radius artery, it is possible 
to know the actual load on the heart by calculating cardiac workload which is obtained from the blood pressure wave- 
form at the proximal portion of the aorta. 

However, while this type of change in the load on the heart may be observed by closely examining the blood pres- 
sure waveform at the proximal portion of the aorta, it is possible to quantitatively express very small changes in the 
waveform by calculatmg the cardiac workload. Accordingly, by obtaining caraiac workload and displaying it along with 
systolic and d.astol.c pressure values, it is possible to carry out an even more precise evaluation of a treatment method 
employing a hypertensive agent. 

The results of calculating cardiac workload for each of the Type I through Type III pulse waveforms are shown in 
rigs. 22 through 24. 

(Modifications) 

The present invention is not limited to the above-described embodiments, but rather a variety of modifications such 
as those below are also possible. For example, a modification may be considered in which the circulatory state param- 
eters are determined without carrying out a measurement of stroke volume SV. 

Namely, from among the circulatory state parameters in this embodiment, inductance L is defined as a fixed value 
and the values of the other parameters are calculated based only on the waveform of the pulsewave measured at the 
test subject's radius artery. In this case, as shown in Fig. 15. it is possible to omit stroke-volume measurer 2 which one 
of the essential elements in the structure shown in Fig. 1. Accordingly, as shown in Fig. 16, when carrying out meas- 
urements in this embodiment, cuff S1 . which is required for the measuring arrangement shown in Fig. 2, is not neces- 
ssry. 

However, when the value of inductance L is defined as a fixed value, the accuracy of the obtained circulatory state 
parameters falls as compared to a method employing the actual measured value of stroke volume. Thus in order to 
compensate for th.s fact, the waveform wi of the radius artery obtained by measurement (i.e.. measured waveform) 
and the waveform W2 of the radius artery obtained by calculation (i.e.. calculated waveform) are superimposed and dis- 
played on output device 6. as shown in Fig. 17. To begin with, the value of inductance L is defined as the aforemen- 
boned fixed value, and a calculated waveform W2 is obtained. This waveform is displayed on output device 6 and the 
degree of coincidence with the measured waveform W1 is observed. Next, the diagnostician determines a suitable 
value which .s drfferent from the fixed value to be inductance L. A calculated waveform W2 is again obtained and the 
degree of coincidence with measured waveform W1 is observed on output device 6. Next, the diagnostician suitably 
determines a number of values for inductance L in the same manner as described above, determines calculated wave- 
forms W2 for each of these inductance L values, and compares each of the calculated waveforms W2 with the meas- 
ured waveform W1 on output device 6. One waveform is selected from among the calculated waveforms W2 which best 
matches measured waveform W1, and the value of the associated inductance L at that time is determined to be the 
optimal value. 

Additionally, note that when modeling the pressure waveform at the proximal portion of the aorta, a squared-off 
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waveform may be employed in place of the triangular waveforms described above. In this case, a waveform results 
which more closely approximately the actual pressure waveform than where approximating using triangular waves. As 
a result, it is possible to calculate even more accurate circulatory state parameters. 

Further, the site of measurement of the pulsewave or stroke volume is not limited to those shown in Figs. 2 and 16. 
5 Rather, any site on the test subject's body may be used. Namely, in the preceding embodiment, cuff S1 was affixed to 
the upper arm of the test subject, however, it may be preferable not to employ a cuff if the convenience of the test sub- 
ject is an issue. 

Therefore, as one example, an arrangement may be considered in which both the radius pulsewave and the stroke 
volume are measured at the wrist. In this type of design, as shown in Fig. 18, a sensor 12 consisting of a sensor for 
10 measuring blood pressure and a sensor for measuring stoke volume are attached to a belt 1 3 of a wristwatch 1 1 . At the 
same time, structural elements 10 other than sensor 12 of a pulsewave analysis device are housed in the main body of 
wristwatch 1 1 . As shown in the figure, sensor 12 is attached in a freely sliding manner to belt 13 via a fastener 14. By 
fastening wristwatch 1 1 to the wrist of the test subject, the sensor presses against the radius artery under a suitable 
pressure. 

is In addition, another embodiment may be considered in which the pulsewave and stroke volume are measured at 

the finger. An example of the structure of a device according to this embodiment in shown in Fig. 19. As shown in this 
figure, a sensor 22 comprising a sensor for measuring blood pressure and a sensor for measuring stroke volume are 
attached to the base of a finger (the index finger in this figure). At the same time, structural elements 10 other than sen- 
sor 22 of the pulse wave analysis device are housed in wristwatch 21 and are attached to sensor 22 via lead wires 

20 23,23. 

Moreover, by combining these two embodiments, other arrangements may be realized such as one in which the 
stroke volume is measured at the wrist and the pulsewave is measured at the finger, or one in which the stroke volume 
is measured at the finger and the radius artery pulsewave is measured at the wrist. 

As in the above-described cases, then, it is possible to realize a structure which does not employ a cuff so that 
25 measurements can be carried out without wrapping an apparatus to the upper arm of the test subject. Thus, the burden 
on the test subject during measurement is reduced. 

On the other hand, a design such as shown in Fig. 20 may be considered which employs a cuff alone. As shown in 
!his figure, a sensor 32 comprising a sensor for measuring blood pressure and a sensor for measuring stroke volume, 
and the structural components 10 other than sensor 32 of a pulsewave analysis device are affixed to the upper arm of 
30 a test subject by means of a cuff. Accordingly, a simpler structure as compared to that in Fig. 2 is realized. 

In addition, in the preceding embodiments, the pulsewave was employed when calculating circulatory state param- 
eters. However, the present invention is not limited thereto; rather other physiological conditions may of course be 
employed for the calculation of the circulatory state parameters. 

35 Claims 



1 . A device for measuring physiological state, comprising: 

a measuring means for measuring physiological state; and 
40 an analysis means for calculating circulatory state parameters including viscoelasticity of an aorta, based on 

the physiological state, as circulatory state parameters which show a circulatory behavior of an arterial system 
from center to periphery of an organism s body. 

2. A device for measuring physiological state according to Claim 1 , characterized in that, additionally, the circulatory 
45 state parameters are resistance in blood vessels at enter of the arterial system due to blood viscosity, inertia of 

blood at the center of the arterial system, resistance in blood vessels at periphery of the arterial system, and vis- 
coelasticity of blood vessels at the periphery of the arterial system. 

3. A device for measuring physiological state according to Claim 2, characterized in that the analysis means is a 
so model having: 

a diode corresponding to an aortic valve; 

a first resistor corresponding to resistance in blood vessels at the center of the arterial system due to blood vis- 
cosity; 

55 an inductance corresponding to blood inertia at the center of the arterial system; 

a first capacitance corresponding to viscoelasticity of the aorta; a second resistor corresponding to resistance 
in blood vessels at the periphery of the arterial system; and 

a second capacitance corresponding to viscoelasticity in blood vessels at the periphery of the arterial system; 
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4. 

w 

5. 

15 



♦ r,K erei H ■ yS ' S meanS calculates the circulatory state parameters by modeling the circulatory 
state of the artenal system using a lumped five parameter model in which a series circuit of the first capaci 
S2L??jn a n connected between a P* ir ° f in P"» ^rminals. a parallel circuit consisting of the second 
,nn S^h!^ ^ °° ^ a ? tenCe iS inS6rted be,ween a P air °t output terminals, and a series circuit consist- 
c^ctenc^ reSiSt ° r iS inS6rted bStWeen 1he 0U1PU, terminal and both termina,s 01 ^ f irst 

A device for measuring physiological state according to Claim 3, characterized in that the physiological state is 
pulsewave at the periphery of the arterial system, and the analysis means determines values of each of elemenS 
making up a lumped f.ve parameter model so that an electric signal corresponding to the waveform of the pulse- 
wave .s obtained from the output terminal when an electric signal corresponding to pressure in the orqanisms left 
cardiac ventricle is provided between the input terminals. organism s left 

A device for measuring physiological state according to Claim 4, characterized in that the electric signal e corre- 
sponding to pressure in the organisms left cardiac ventricle may be approximately expressed by the sinusoidal 

W3V6. 
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e = E m 'sin (tt/1 s )t 

where E m ' is systolic pressure, t s is the time duration of rising pressure in the left cardiac ventricle, and t is time. 

A device for measuring physiological state according to Claim 1 , characterized in the provision of a blood pressure 
Sra U m^e?s meanS Ca,CU,atin9 the preSSUre waveform at the organism's aorta based on the circulatory state 

A device for measuring physiological state according to Claim 6, characterized in that the circulatory state param- 
th™? re f;J anC f ,n f s blood ves5els ai ihe cenier of the arterial system due to blood viscosity, blood inertia at 
tetnrZ Z SL w *1 T^'. resistance in the blood vessels at the periphery of the arterial system, and viscoe- 
lasticity of the blood vessels at the periphery of the arterial system. 

^^XT^IX* 09 ™ 51316 according to Claim 7l characteri2ed in that the Wood 

a diode corresponding to the aortic valve; 

a first resistor corresponding to resistance in blood vessels at the center of the arterial system due to blood vis- 
cosiiy, 

an inductance corresponding to blood inertia at the center of the arterial system; 
a first capacitance corresponding to viscoelasticity of the aorta; 

a second resistor conesponding to resistance in blood vessels at the periphery of the arterial system and 
a second capacitance corresponding to viscoelasticity in blood vessels at the periphery of the arterial system- 
wherein the blood pressure calculating means determines the circulatory state parameters by modeling 
the circulatory state of the arterial system using a lumped five parameter model in which a series circuit of the 
irst capacitance and the diode is connected between a pair of input terminals, a parallel circuit consisting of 
the second capactence and the second resistor is inserted between a pair of output terminals, and a series 
crcurt consis ting of the inductance and the first resistor in inserted between the output terminal and both ter- 
minals of the first capacitance, and defines waveform of voltage between both terminals of the first capacitance 
as pressure waveform at the aorta. ^ 

9. A device for measuring physiological state according to Claim 8, characterized in that the physiological state is the 
pulsewave at the periphery of the arterial system, and the blood pressure measuring means determines the values 
of each of the elements making up the lumped five parameter model so that an electric signal conesponding to the 
waveform of the pulsewave is obtained from the output terminal when an electric signal corresponding to pressure 
m the organism s left cardiac ventricle is provided between the input terminals. 

10. A device for measuring physiological state according to Claim 6 or Claim 7. characterized in that the physiological 
ISS™ hT S6wava at peri 0 her * 01 the arte " al astern, the device has a distortion calculating means for cal- 
cu atmg distortion in the pulsewave from the waveform of the pulsewave. and the blood pressure measuring means 
determines the circulatory state parameters based on the correlation between the circulatory state parameters and 
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10 



the distortion in the pulsewave. 

11. A device for measuring physiological state according to one of Claims 6 through 9, characterized in the provision 
of a stroke-volume detector which detects the stroke volume in the organism, and in that the blood pressure calcu- 
lating means adjusts the value of the circulatory state parameters so that the calculated value of stroke volume 
obtained from the pressure waveform at the aorta and the actual measured value of stroke volume measured by 
the stroke-volume measurer are equivalent. 

12. A device for measuring physiological state according to one of Claims 6 through 9, characterized in the provision 
of a workload calculating means for calculating workload on an organism's heart based on the pressure waveform 
at the aorta 

13. A device for measuring physiological state, comprising: 

15 a measuring means for measuring physiological state at the periphery of an arterial system; 

a blood pressure calculating means for determining circulatory state parameters, which display the circulatory 
state of the arterial system, from physiological state, and calculating pressure waveform at an aorta based on 
the circulatory state parameters; and 

a ventricle state calculating means for calculating the ventricular contraction phase, ventricular expansion 
20 phase, and duration of rising pressure in a ventricle. 

14. A device for measuring physiological state according to Claim 6, characterized in that the measuring means meas- 
ures the pressure waveform of the pulse at the periphery of the body, and further characterized in the provision of 
a comparing means for comparing the pressure waveform of the pulse at the periphery of the body and the pres- 

25 sure waveform at the aorta, and a diagnosing means for carrying out a diagnosis based on results of the compari- 
son. 

15. A device for measuring physiological state according to Claim 14, characterized in that the measuring means 
measures a waveform at a radius artery of the body, the comparing means determines the difference in maximum 

30 value, minimum value and average value between the waveform of the radius artery and the pressure waveform at 
the aorta, and the diagnosing means carries out a warning when the difference obtained is within a specified range. 
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FIG. 9 



(I CALCULATI ON OF n Alan m h 
* 



eps<— 100 

I 



a ' *-3 
I 



✓ S201 
S202 




NO 



NO 



|l CALCULAT|ON OF cu ~~p ^ 

< ( Smax- I V p (t/)- yi | 

+YES ~ 

S nax"- I Vp(t/)- yi | 

< ( s m ax"- I VpWj-yg j 
f +YES 

j S max*~ I Vp(t 2 Q-y 2 j 



S203 
S204 





S205 
S206 




<C^fnax<ePSl 


iYES 












S max*~eps1 




* 




Qf ' «— a ' +0. 1 




S212 
S213 
YES 



42 



EP 0 818 175 A1 



FIG. 10 
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FIG. 11 
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FIG. 18 
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FIG. 22 
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FIG. 25 
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FIG. 27 
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FIG. 32 
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